Evidence for universal structure in galactic halos 
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The late infall of dark matter onto a galaxy produces structure (such as caustics) in the distribution 
of dark matter in the halo. We argue that such structure is likely to occur generically on length 
scales proportional to I ~ toVrot, where to is the age of the universe and Urot is the rotation velocity 
of the galaxy. A set of 32 extended galactic rotation curves is analyzed. For each curve, the radial 
coordinate is rescaled according to r — > f = r(wo/wrot), where we choose no = 220km/s. A linear fit 
to each rescaled rotation curve is subtracted, and the residuals are binned and averaged. The sample 
shows significant features near f = 40 kpc and f = 20 kpc. This is consistent with the predictions of 
the self-similar caustic ring model of galactic halos. 
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The rotation curves of most spiral galaxies are approx- 
imately flat, i.e. {v{r) ~ constant) for r much larger than 
the disk radius [y, where v{r) is the circular velocity of 
gas at radial coordinate r. The flatness of galactic rota- 
tion curves implies that galaxies are surrounded by halos 
of dark matter, and that the halo density falls off ap- 
proximately as pDMif) oc 1/r^. The curves are far from 
exactly flat, however. They have all sorts of irregulari- 
ties which may be referred to as "bumps" . When a bump 
occurs at r >> disk radius, it indicates structure in the 
dark matter distribution, i.e. a deviation from a perfect 
Pdm oc l/r^ law. We are motivated by the possibility 
that such structures have universality and that a statis- 
tical analysis of well measured extended rotation curves 
may reveal that. 

In refs. [piP], 32 extended galactic rotation curves were 
selected and analyzed to test the validity of a modifi- 
cation of Newtonian dynamics (MOND) as an explana- 
tion of the "dark matter problem" , i.e. the discrepancy 
between the flatness of galactic rotation curves and the 
exponential fall-off of galactic luminous matter distribu- 
tions. We analyze the same set of rotation curves with a 
different goal, namely to search for structure within dark 
matter halos. We adopt the conventional wisdom that a 
rotation curve is a tracer of the galactic mass distribution 
according to the laws of Newtonian dynamics. However, 
the selection criteria used by the authors of refs. gj^ - 
that each rotation curve is an accurate tracer of the ra- 
dial force law, and that it extends far beyond the edge 
of the luminous disk - are appropriate from the point of 
view of our analysis. We do not know of another data 
set of comparable quality. We use the complete data set 
described there, without any cuts of our own. 

One broad argument why structure may be present in 
galactic halos is that the halos form as a result of the 
infall of the dark matter surrounding the galaxy [H. If 
the dark matter particles are coUisionless, they oscillate 
back and forth numerous times before they are virialized 
by inhomogeneities in the galactic mass distribution [H. 



The non- virialized flows of dark matter produce structure 
within the halo, e.g. the caustics described below. 

If structure is present, one would expect it to occur 
on length scales set by t Vrot where t is the age of the 
galaxy and Vrot is its rotation velocity. Indeed Vrot sets 
the scale for the velocities with which all dark matter 
particles move about in the halo. Since all galaxies have 
nearly the same age, we expect that any regularity shared 
by the different galaxies of a set will reveal itself most 
readily after all lengths have been rescaled according to 

We are motivated in large part by the fact that caus- 
tics form in the non-virialized flows associated with the 
late infall of dark matter onto a galaxy. By definition, a 
caustic is a place in physical space where the dark matter 
density is very large. In the limit where the infalling dark 
matter has zero velocity dispersion, the dark matter den- 
sity diverges at the caustics. Two types of caustic form: 
mner and outer [p|. The outer caustics are 2-spheres sur- 
rounding the galaxy. One such caustic is located wher- 
ever an outfiow of dark matter reaches its maximum ra- 
dius before falling back in. The inner caustics are rings. 
One caustic ring is located wherever the particles with 
the most angular momentum in a given infiow reach their 
distance of closest approach to the galactic center before 
moving back out. 

Outer caustic spheres and inner caustic rings are ex- 
amples of the sort of structure we are looking for. Caus- 
tic rings are the more likely candidate for detection in 
an analysis of galactic rotation curves because they have 
the higher density contrast and because they are located 
closer to the galactic center, where the rotation curves 
are measured. Furthermore, a specific proposal [Q has 
been made for the radii a„ of caustic rings: 
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where h is the present Hubble rate in units of 



100km/(s Mpc), and jmax is the maximum of the dimen- 
sionless angular momentum distribution of the infaUing 
dark matter particles in the self-similar infall model of 
galactic halo formation. 

The infall of dark matter is called self-similar B if it 
is time-independent after all distances are rescaled by a 
time-dependent scale R{t) and all masses are rescaled by 
the mass M{t) interior to R{t). Usually R{t) is taken to 
be the turn-around radius, defined as follows. Consider 
all particles which are about to fall onto the galaxy for 
the first time in their history at time t. Such particles 
are said to be at their 'first turn-around'. In a spher- 
ical model, they have zero radial velocity then. Their 
distance to the galactic center is the turn-around radius 
R{t) at that time. In the case of zero angular momen- 
tum and spherical symmetry, the infall is self-similar if 
the initial overdensity profile has the form 
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where Mq and e are parameters [Bl . e must be in the range 
< e < 1. The rotation curve is flat if < e < 2/3 |||. In 
models of large scale structure formation, e is predicted 
to lie in the range 0.2 to 0.35 M. The self-similar infall 
model was generalized in ref. |^ to include the effect of 
angular momentum. It was found that, in addition to the 
condition on the initial overdensity profile, self-similarity 
requires the angular momentum distribution i(t) to have 
the time-dependence 
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where j is a dimensionless and time-independent distri- 
bution. The maximum of the j distribution is the quan- 
tity jmax which appears in Eq. (1). Good agreement of 
the self-similar model with the properties of our own 
galaxy was found H for parameter values e = 0.20 to 
0.35, j ~ 0.2 and h ~ 0.7 where j is the average of the 
j distribution. 

Eq. (1) assumes e = 0.3. For e = 0.2, 



oi ~ 36kpc 
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but the ratios a„/ai are almost the same as in the e = 0.3 
case. The ratios happen to be close to a„/ai = 1/n over 
the range (0.2 < e < 0.35) of interest. 

Eq. (1) predicts the caustic ring radii of a galaxy in 
terms of its first ring radius ai. If the caustic rings lie 
close to the galactic plane they cause bumps in the rota- 
tion curve at the caustic ring radii. As a possible example 
of this effect, consider the rotation curve of NGC3198, 
one of the best measured and a member of the set selected 
in refs. HJ^. It has three faint bumps at radii: 28, 13.5 
and 9 kpc, assuming h = 0.75. The ratios happen to be 



consistent with Eq.(l) assuming the bumps are caused 
by the first three (n = 1, 2, 3) ring caustics of NGC3198. 
Moreover, since Uiot = 150 km/s, j,„ax is determined in 
terms of e. For e = 0.3, jmax = 0.28. Note that the 
uncertainty in /i is a systematic effect that can be cor- 
rected for when determining jmax because the bump radii 
scale like 1/h' where h' is the Hubble rate assumed by 
the observer in constructing the rotation curve, and the 
caustic ring radii scale as 1/h. Rises in the inner rotation 
curve of the Milky Way were also interpreted m as due 
to caustics n = 6,7,8,9,10,11,12 and 13. This deter- 
mined the value of jmax of our own galaxy to be 0.263 
for e = 0.3, the value we assume henceforth. The first 
five caustic ring radii in our galaxy are then predicted to 
be: 41, 20, 13.3, 10, 8 kpc. 

According to the self-similar caustic ring model, each 
galaxy has its own value of jmax- Over the set of 32 
galaxies selected in refs. [||j^, jmax has some unknown 
distribution. However, the fact that the values of jmax of 
NGC3198 and of the Milky Way happen to be close to one 
another, within 7%, suggests that the jmax distribution 
may be peaked near a value of 0.27 . Our strategy is to 
rescale each rotation curve according to 
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and to add them together in a way made precise be- 
low. Since Eq. (1) predicts the nth caustic radius a„ to 
be distributed like jmax for all n, and it fixes the ratios 
ttn/ai ~ 1/n, the sum of rotation curves should show the 
jmax distribution, once for n = 1, then at about half the 
n = 1 radii for n = 2, then at about 1/3 the n — I radii 
for n = 3, and so on. If the jmax distribution is broad, 
the sum of rotation curves is unlikely to show any fea- 
ture. However, if it is peaked, then the sum should show 
a peak for n = 1 at some radius, then again at 1/2 that 
radius for n = 2, at 1/3 the radius for n = 3, and so on. 
If the jmax distribution is peaked at 0.263 (the value for 
the Milky Way) , the peaks in the sum of rotation curves 
should appear at 41 kpc, 20 kpc, 13.3 kpc, and so forth. 
We search for overdensities in the halo dark matter 
distribution as manifested by upward fluctuations in ve- 
locity relative to the underlying flat rotation curve of the 
galaxy. Unfortunately, rotation curves are never exactly 
fiat. They often rise or fall systematically with increas- 
ing radius, and all show a sharp drop in rotation velocity 
near the galactic center. Furthermore, individual sources 
(i.e. gas) in a galaxy have peculiar velocities indepen- 
dent of any caustic structure in the galaxy. We require 
a definite procedure for subtracting out the background 
rotation of the galaxy, and for quantifying the "noise" 
due to the peculiar velocities of sources. The procedure 
we adopt is as follows. Each rotation curve is expressed 
as a set of radii rt in kpc and rotation velocities Vi in 
km/sec. Each (r^, Vi) is given by a data point in the 



rotation curves published in refs. ||,^-[| Now, in the case 
of a flat rotation curve, Urot is just the average velocity 
of the outer rotation curve, which we denote v. In the 
case of a rising or falling rotation curve, the "overall" 
rotation velocity Wrot is less precisely defined. We adopt 
the prescription Wrot = v in all cases. Thus the rescaled 
radii f,- are 
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where v is calculated according to the following proce- 
dure. We are interested in the outer, approximately flat 
portion of the rotation curve, so we remove the inner 
points by specifying a cut: all points with rescaled radii 
fi < lOkpc are thrown out. The remaining points in 
the rotation curve are then fitted to a line, and v is the 
corresponding average. A subtlety arises because the av- 
erage velocity v is calculated after the cut is applied, but 
the cutoff is defined in rescaled coordinates and therefore 
depends on v. The cutoff procedure is in fact performed 
iteratively: we start with a guess for v, calculate the cut- 
off radius, recalculate v based on the cutoff radius, and 
so on until both quantities converge (usually after one or 
two iterations). 

Fitting the outer rotation curve to a line may seem an 
arbitrary choice. In practice, it works quite well: devi- 
ations from a linear fit are typically less than lOkm/s 
in galaxies with a typical rotation velocity of 200 km/s. 
However, to verify the robustness of our conclusions, we 
also performed the analysis with quadratic polynomial 
fits to the rotation curves, with no substantial change in 
the results. (Of course, a fit to a high enough order poly- 
nomial will remove any features present in the rotation 
curve!). 

Once the linear fit is calculated, this background ro- 
tation is subtracted off the rotation curve, leaving a set 
of peculiar velocities Svi. It is then straightforward to 



calculate an rms "noise" y/{5v 2) for each galaxy. 
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where N is the number of points and Nm is the number of 
degrees of freedom in the fit {Nfn = 2 for a fit to a line). 
We adopt y^{5v ^) as the size of the error bar on the 6vi 



* These curves are heavily processed and do not represent the 
velocities of individual sources. See, for example, Ref. ||l(l] 
and references therein for a discussion of the observational 
details relating to one galaxy in our sample, IC2574. In the 
case of the galaxy M33, whose published rotation curve |Q,m| 
is so densely sampled that the data points overlap with one 
another, we choose to sample the published curve at intervals 
of 0.25 kpc, a sufficient resolution for our purpose here. 



for a given galaxy. This error is a measure of the intrinsic 
peculiar velocities of the sources, which we assume to be 
random. It is considerably more conservative than the 
quoted observational errors on the points in the rotation 
curve. Finally, the peculiar velocities are expressed in 
dimensionless units Svi = Vi/y^{Sv ^), and the sample of 
galaxies is averaged into radial bins: 
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where Ni is the number of data points in the z-th bin. The 
assigned error on each bi is then simply 1/ y/Wl. Figure H 
shows the complete set of 32 galaxies averaged into 2 kpc 
bins. 

There are two features evident at roughly 20 and 
40 kpc. A fit to two Gaussians (with amplitude, width 
and mean left as free parameters) plus a constant indi- 
cates features at 19.4 ± 0.7 kpc and 41.3 ± 0.8 kpc, with 
overall significance of 2.4cr and 2.6cr, respectively. Fig- 
ure also shows the fitted curve. When the same fit is 
applied to the same data in 1 kpc bins, the significance 
of the two peaks is 2.6cr and 3. Oct respectively. The lo- 
cations of the features agrees with the predictions of the 
self-similar caustic ring model with the jmax distribution 
peaked at 0.27. The use of Gaussians to fit the peaks in 
the combined rotation curve was an arbitrary choice in 
the absence of information on the jmax distribution. 
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FIG. 1. Binned data for 32 galaxy sample, with peaks fit 
Gaussians 

We have shown that there is evidence for universal 



structure in the dark halos of spiral galaxies, in that 
the rotation velocity has a statistically significant ten- 
dency to fluctuate upward at rescaled radii f ^ 20 and 
40 kpc. While the significance of these peaks in the data 
set analyzed is not overwhelming, roughly 2.5a for each 
peak when considered with appropriately conservative er- 
ror bars, our analysis provides a tantalizing indication 
that the structure is real and that it is within observa- 
tional reach of a more detailed survey of galactic rotation 
curves. 

It is particularly striking that the positions of the peaks 
coincide with the radii of the n — I and n — 2 caustic 
rings in the self-similar infall model for jmax ~ 0.27. This 
strongly suggests that the jmax distribution is indeed 
peaked near the value of 0.27, as the earlier estimates of 
the j„rax values of NGC3198 (0.28) and the Milky Way 
(0.265) had suggested. 
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